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ABSTRACT 



Context. Owing to the limited spatial resolution and the weak polarization signal coming from the quietest regions on 
the Sun, the organization of the magnetic field on the smallest scales is largely unknown. 

Aims. We obtain information about the magnetic flux present in the quiet Sun by comparing radiative MHD simulations 
with observations, with particular emphasis on the role of surface dynamo action. 

Methods. We synthesize Stokes profiles on the basis of the MHD simulation results. The profiles are degraded taking 
into account the properties of the spectropolarimeter (SP) on board of the Hinode satellite. We use simulation runs 
with different magnetic Reynolds numbers (Rm) and observations at different heliocentric angles with different levels 
of noise. 

Results. Simulations with an imposed mixed-polarity field and Rm below the threshold for dynamo action reproduce the 
observed vertical flux density, but do not display a sufficiently high horizontal flux density. Surface dynamo simulations 
at the highest Rm feasible at the moment yield a ratio of the horizontal and vertical flux density consistent with 
observational results, but the overall amplitudes are too low. Based on the properties of the local dynamo simulations, 
a tentative scaling of the magnetic fleld strength by a factor 2 — 3 reproduces the signal observed in the internetwork 
regions. 

Conclusions. We flnd an agreement with observations at different heliocentric angles. The mean fleld strength in 
internetwork, implied by our analysis, is roughly 170 G at the optical depth unity. Our study shows that surface dynamo 
could be responsible for most of the magnetic flux in the quiet Sun outside the network given that the extrapolation 
to higher Rm is valid. 

Key words. Sun: granulation. Sun: photosphere 
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1. Introduction 

The origin of the sma ll scale magnetic flu x found in the 
quiet Sun is uncertain (|de Wiin et al.|[2Q08[ ). The estimated 
« order of magnitude for the magnetic Reynolds number of 
I the granulation flow indicates that a substantial part of 
. the magnetic field in the quiet Sun could be generated lo- 

■ cally through dynamo action driven by near- surface con- 
[ vective flows (jPetrovav fc Szakalvl 1 19931 ICattaneol 1 19991 
' IVogler fc Schiisslerll200^ The simulations suggest that the 

■ magnetic field could be organized in mixed-polarity struc- 
tures down to very small spatial scales. The simulations 
exhibit a mostly horizontal field in the photospheric layers 
in the for m of low-lying loops conn ecting nearby opposite 
polarities (|Schiissler fc Voderl [20081 ). 

Substantial observational evidence has been gathered 
that internetwork flux is dominated by strongly inclined, 
almost horizontal magne tic fields ([Orozco Suarez et al.l 
l2007allbl : iLites et al.l I2008D . The v alidity of the d e duced 
properties has been questioned by lAsensio RamosI (|2009[ ) 
who argues that the influence of noise has not been ad- 
equately taken into account in such studies. Nevertheless, 
several author s report loop- l ike horizontal field structures of 
different sizes (|Martinlll988[lMartmez Gonzalez et al.ll2007l : 



Harvev et all l2007l : iGenteno etld] 120071 : llshikawa et"all 
20081) . Such structures could be due to local dynamo ac- 



tion, but small-scale flux emergence (|Cheung et al.ll2008h 
or flux expulsion of a pre-existing field by granular flows 
(|Steiner et al.|[2QQ8l ) probably also contribute to the hori- 
zontal fiux. 

Another open question concerns the amount of mag- 
netic flux contained in the internetwork. The observa- 
tions give a wide range of values that vary with the 
spatial resolution and the diagn o stic technique used (e.g , 
iDommguez Cerdefia et all (I2003D: iKho menko et al.' f20^; 
feerdyugina fc Fluril (|2Qq1 ): lAsensio Ramos et al . (2 0q3)) . 
From many results obtained (see Ide Wiin eFal.l 120081 
ISolankil l2009l for overviews), we mention here only those 
that directly concern the work presented in this paper. 
Comparing the center-to-limb observations of the scatter- 
ing polarization in the Sr I 4607 A line with the sig- 
nal synthesized from the 3 D hydrodynamical simulations, 
iTruiillo Bueno et al.l (|2004{ ) inferred that the mean strength 
of the internetwork field is (B) ~ 100 G, under the as- 
sumption that the magnetic field is isotropically tangled 
at subresolution scales and that i t fills the whole res- 
olution element. iKhomenko et al.l (|2005al ), on the other 
hand, compared the observed Stokes V amplitudes of the 
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Fig. 1. Maps of longitudinal {B^^^) and transversal [B^^^) apparent flux density for the observed data sets I (left) and 
III (right) from Table [H The small squares indicate the sizes of the MHD snapshots. They also mark the position of the 
regions to which synthesized maps from simulation data are compared (Figures H] and [5j). 



visible and infrared Fe I lines with the profiles synthe- 
sized from 3D radiative MH D simulations, and co ncluded 
that {B) ^ 20 G (see also iBello Gonzalez et all (|2009)). 
I Sanchez Almeida et al.l (|2QQ3h used Boussinesq 3D simula- 
tions of local dynamo action to suggest that a snapshot with 
mean longitudinal magnetic field of 50 G can reproduce 
both the observed Hanle and Zeeman signals. However, the 
simulations used by these authors are rather idealized and 
the solar atmospheric structure was arbitrarily introduced 
for comparison with the observational data. 

In this paper, we use 3D radiative MHD simulations of 
the solar photosphere to obtain an estimate of the true 
magnetic fiux densi ty in the quiet S un. O ur approach 
differs from that of iKhomenko et al.l (|2QQ5al ) in that we 
also consider fields produced by a local dynamo action. 
Also, the synthesized polarization signals are compared 
with the highe r resolution data , obtained with the spec- 
tropolarimeter (|Lites et al.ll2QQlf ) of the Solar Optical tele- 



scope (|Tsuneta et al.ll2QQ8l ) on board of the Hinode satellite 
(jKosud et al.ll2QQ7D . The comparison with the simulations 
at the highest feasible at the moment also allows us to 
gain insight into the role of the local dynamo action for the 
quiet Sun magnetism. 

2. Observations 

We consider three data sets obtained with the spectropo- 
larimeter (SP) on board Hinode. Details are given in 
Table [H Data sets I and II were recorded near disc center, 
while data set HI was obtained closer to the limb. Data sets 
I and III were obtained in the scan mode of the Hinode/SP, 
with an exposure time of 9.6 s per slit position. They cover 
quiet Sun regions at the disk center and a region near the 
south solar pole, respectively. In the case of data set HI, 
we use a 32'^ wide strip perpendicular to the scan direc- 
tion that corresponds to /i = 0.4 (/i being the cosine of the 
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Table 1. Details of the three Hinode/SP data sets used. The last three columns, from left to right, give the mean values 
of the transversal (with respect to the line of sight), unsigned longitudinal and signed longitudinal magnetic flux density, 
respectively. 





date 


time 


FOCUS 




size 


exposure 


noise 


i-Bapp) 


{\Bapp\) 


i-Bapp) 




[dd/mm/yy] 


UT 






[arcsec] 


time [s] 


10-^/c 


Mx/cm^ 


Mx/cm^ 


Mx/cm^ 


set I 


10/09/07 


08:00:00 


2031 


1 


58 X 164 


9.6 


0.8 


56.7 


9.7 


1.6 


set II 


27/02/07 


00:20:09 


2048 


1 


0.16 X 164 


67.2 


0.25 


54.6 


11.0 


1.7 


set III 


09/09/07 


13:05:05 


2029 


0.4 


160 X 32 


9.6 


1 


53.5 


6.8 


0.3 



Table 2. The simulation parameters. Given are the size 
of computational domain and vertical (Sz) and horizontal 
{6x) grid spacing. 



Run 


height/width 


6z/6x 


Rm 




[Mm] 


[km] 




mixed polarity 


1.4/6.0 


14/20.8 


- 300 


dynamo C 


1.4/4.86 


10/7.5 


- 2600 


dynamo G 


1.4/4.86 


7/5 


- 5200 



heliocentric angle). Corrections for various instrumental ef- 
fects are made using the procedure sp.prep, included in 
the Solar SoflQ package. The procedure gives th e longitudi- 
nal a nd transversal magnetic flux density maps (|Lites et al.l 
I2008D shown in FiglH The mean values are given in Table [H 
The rms continuum contrast values for data sets I and III 
are 7.5% and 5.1% respectively. 

Data set II is a time series recorded with fixed slit po- 
sition. It consists of 103 scans at solar disk center, each 
with an exposure time of 9.6 s. After applying a temporal 
running mean, the effective exposure time becomes 67.2 s, 
which gives a significantly lower l evel of noise . This data 
set has previously been us ed by iLites et aD (|2008[ ) and 
lOrozco Suarez et all (|2007b[ ). 

3. Simulation data 

We use results from 3D radiative MHD simulations of a 
layer containin g the s olar s urface, carrie d out with the 
MURAM code (|VQder|[2QQl IVoder et al.ll2005[ ). Non-grey 
LTE radiative transfer and partial ionization are taken into 
account. We compare snapshots from several simulation 
runs. The basic properties of the computational domains 
are given in Table [2j In all runs, the top of the simulation 
box is located about 500 km above the average height level 
of optical depth unity. The side boundaries are periodic, 
whereas the bottom boundary is open, permitting free in 
and outflow of matter. The magnetic field is vertical at the 
top and bottom boundaries. 

The first run, henceforth referred to as the 'mixed po- 
larity' run, simulates the decay of the magnetic field in a 
mixed polarity region. In this run, local dynamo action 
(|VQgler fc SchiissTerl I2OO7I ) does not occur since the mag- 
netic Reynolds number is below the threshold for dynamo 
action. The run starts with a vertical magnetic field of 
\B\ = 200 G, in a checkerboard-like 2x2 pattern, with 
opposite polarities in adjacent parts. The field is concen- 
trated and redistributed by the convective motions; the op- 
posite polarities are pushed together, which results in flux 
cancellation and almost exponential decrease of the mean 
magnetic field strength. The snapshots with {\B\) = 35 G 

^ http:/ /www. lmsal.com/solarsoft/ 



and {\B\) = 20 G averaged over the surface r = 0.1 
are used in our s tudy. T he same snapshots were used by 
iKhomenko et al.l (|2QQ5al lbh. 

The second and third group of snapshots are taken from 
the runs with a magnetic Reynolds number of the flow suf- 
ficiently high for small-scale turbulent dynamo action to 
take place. A seed field of |5o| = 10 mG grows exponen- 
tially in time until the sat uration level is reached. In the 
dynamo run C, described in lVogler fc Schiisslerl (|2QQ7f ), the 
mean magnetic field strength at this phase is = 23 G 
at r = 0.1. In the case of the dynamo run G, the higher 
magnetic Reynolds number leads to a saturation level of 
the magnetic energy which is a factor of 1.7 higher than 
in run C. Figure [2] shows the mean magnetic field strength 
(averaged over surfaces of constant r) as a function of the 
optical depth for one snapshot from each run. The field 
from dynamo C run is multiplied by a factor 1.5, in order 
to illustrate the similarity of the scaled optical depth pro- 
files in the dynamo runs. This factor differs only by ~ 15% 
from \/T77, the square root of the corresponding total mag- 
netic energy ratio. This indicates that a tentative scaling 
of the dynamo-generated field to represent higher Reynolds 
numbers is not completely unreasonable. The dashed lines 
show the optical depth profiles of the horizontal magnetic 

field, i.e., (Bhor) = {\J Bl + B'^) and the dotted lines rep- 
resent the average vertical field {\Byert\)- The dynamo runs 
show a significantly larger {Bhor) than {\Bvert\)- Their ratio 
reaches values between 2 and 4 in the optical depth inter- 
val —2 < logr < — 1. The mixed polarity snapshot, on the 
other hand, has a (Bhor) and {\Byert\) of similar magnitude 
over all heights. 

The last group of snapshots is taken from a simulation 
run that continues dynamo run C with a superposed unipo- 
lar vertical field (green lines in Figure [2]). Such a superpo- 
sition might be a way of describing a weak network region. 
The mean vertical field strength is thereby increased to 
around 36 G at r = 0.1 and the strength of the horizontal 
field to around 50 G at r = 0.1. 



4. Spectral synthesis 

The simulation results have been us ed as input for the 
LTE radiative transfer code SPINOR (|Frutiger et al.ll2QQQ[ ) 
to synthesize the Stokes profiles for the heliocentric an- 
gles ^ = 0° and 66° (corresponding to /i = 1 and 
fi = 0.4, respectively). The spectral range that contains 
the Fe I lines at 630.15 and 630.25 nm is sampled in steps 
of 7.5 mA. The Fe abundanc e used in the synthesis has been 
taken from lTheveninI (|l989l ) and the values of the oscillato r 
strengths from the VALD database (|Piskunov et al.|[T995f ). 
We then applied a re alistic point spread function (PSF, 
iDanilovic et al.l (|2QQ8f )) to the maps of synthesized Stokes 
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Fig. 3. Histograms of the mean horizontal (left), the unsigned vertical (middle) and the signed vertical (right) magnetic 
flux density calculated for 1" x 7" size regions covered by dataset I. The vertical lines mark the mean values for the whole 
field of view. 
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Fig. 2. Average total (solid), horizontal (dashed) and ver- 
tical (dotted) mean magnetic field strength from MURaM 
simulations as a function of the optical depth. Snapshots 
from different runs are indicated by line color: mixed po- 
larity (blue), dynamo G (red), dynamo C multiplied by fac- 
tor of 1.5 (black) and dynamo C with background unipolar 
field (green). 

profiles. The PSF takes into account the basic optical prop- 
erties of the Hinode SOT/SP system and a small defocus. 
Applying an ideal PSF without defocus reduces the orig- 
inal rms contrast values of the simulated continuum map 
from 14.5% to 8.5% {0 = 0°) and 11.2% to 5.5% {0 = 66°), 
respectively. For the comparison with the observations at 
6> = 0°, we used a defocus of —1.5 mm, which degrades the 
continuum contrast of the simulation to the observed value 
of 7.5%. A value of —0.75 mm is used for the synthesized 
data at 6> = 66° in order to match the continuum contrast of 
5.1% deduced from dataset III. The difference in the focus 
of the SOT between data sets I and III amounts to approx- 
imately 2 steps of the SOT focus mechanism. Taking into 
account that there is an uncertainty of one focus step in 
the focusing mechanism, the amounts of defocus obtained 
from reproducing the observed data sets seem plausible. 

The appropriate PSF is applied to the 2D maps at each 
wavelength position for every Stokes parameter. The de- 
graded maps are then rebinned to the pixel size of the 



Hinode/SP detector. To take into account the spectral res- 
olution of the spectropolarimeter, the profiles are convolved 
with a Gaussian function of 25 mA FWHM and resampled 
to a wavelength spacing of 21.5 mA. Next, a noise level 
correspondin g to the observatio ns is added, and finally the 
procedure bv lLites et all (|2QQ8f ) is used to calculate the lon- 
gitudinal and transversal apparent magnetic flux densities, 
^app and B^^^. 

5. Results 

We present our results in the following sequence: (1) snap- 
shots from the simulation runs are compared with the disc 
center observations in terms of the mean values of and 

app 

B^pp] (2) the spatial distribution of the magnetic features in 
the observed and in the synthesized maps is compared; (3) 
a comparison of the probability density functions (PDFs) 
from the simulations and from the low-noise data set II is 
presented. Finally, we show how the mean magnetic flux 
density changes with the heliocentric angle and compare 
the synthesized and observed signal at a heliocentric angle 
of l9 = 66°. 

5.1. Comparison with the observations at disc center 

When comparing the observational signals with the ones 
synthesized from the simulations, one has to take into ac- 
count that the solar surface area covered by the observa- 
tions of data set I is much larger than that comprised by 
the simulations. The small squares in Fig [1] indicate the 
actual size of the simulation snapshots. In order to take 
into account this proportion, we divide the region covered 
by the observations into subdomains of 7" x 7" size. For 
each of these we calculate the mean transversal (B^pp), the 
mean unsigned longitudinal {\B^pp\) and the mean signed 
longitudinal (B^pp) apparent magnetic flux densities. The 
histograms of these quantities are shown in Fig [3l By con- 
sidering the contribution of each subdomain separately, we 
can exclude the contribution of the network, which extends 
the wings of the histograms towards higher values of {B^pp) 
and (l^^^^l). The maxima of the distributions, on the other 
hand, give an estimate of the typical magnetic flux den- 
sity value in the internetwork regions. Consequently, when 
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Table 3. Mean apparent magnetic flux densities from observations and simulations. The values in parentheses are the 
mean values at the original resolution of the simulations, before spatial smearing and addition of noise. 



Run 




\\-^app\) 


/ R-^ \ 
\^app) 




\J^app)/[\-Dapp\) 


observations 


55 


a 
D 


o 
-z 




n o 
9.Z 


mixed polarity (30 G) 


39 (24) 


13 (19) 


-0.2 (-0.3) 




3.0 (1.2) 


mixed polarity (20 G) 


37 (15) 


6.5 (9.0) 


-0.6 (-0.3) 




5.7 (1.8) 


dynamo C 


36-37 (21-28) 


2.7-3.0 (6.4-7.8) 


-0.1-0.0 (-0.2- 


0.0) 


12-13 (3.3-3.6) 


dynamo C (H-uni. field) 


45-51 (51-65) 


9-15 (18-27) 


9-14 (9-16) 




3.4-5 (2.4-2.8) 


dynamo G 


39 (39) 


3.0 (9.6) 


0.1 (0.1) 




13 (4.1) 


noise 


36 


2.3 


0.0 






dynamo C (mf — 3) 


51-60 (70-81) 


5.3-6.0 (19-24) 


-0.1-0.0 (-0.1 - 


0.0) 


9.6-10 (3.4-3.7) 


dynamo G {mf — 2) 


52 (77) 


4.7 (19) 


0.2 (0.2) 




11 (4.0) 



the whole observed region is taken into account, the mean 
values (marked by the vertical lines in Fig [3]) are higher 
than the values retrieved from the maxima. The mean val- 
ues over the wh ole FOV are in agreement with the val- 
ues obtained by iLites et all (|2QQ8i) and show a ratio of 

The values retrieved from the distribution maxima are 
given in the first row of Table [3l The other numbers are the 
results obtained from the synthesized Stokes profiles from 
the simulation snapshots; values determined from the snap- 
shots at their original resolution are given in parenthesis. 
They roughly correspond to the mean vertical and horizon- 
tal field strengths near r = 0.1 (cf. Figure [2]). Thus, the 
ratio (B^pp) / (B^pp) obtained from the simulations at the 
original resolution reflects the ratio of the underlying mag- 
netic fields at r = 0.1. The values are close to unity for the 
mixed-polarity run and about 3 for the dynamo snapshots. 
For the dynamo snapshots, the values {\B^pp\) after spatial 
smearing and application of noise are considerably reduced 
compared to the noise- free, unsmeared case owing to the 
presence of mixed polarities on very small scales. On the 
other hand, (B^pp) is increased for all groups of snapshots, 
except for the dynamo run C with unipolar background 
field. This is the result of the noise which is added to the 
Stokes profiles in order to simulate the Hinode/SP obser- 
vations. The row labeled "noise" in Table [3] gives values 
determined from pure white noise with a standard devia- 
tion corresponding to the noise level of data set I. We give 
the mean values from 100 realizations. The value of (B^pp) 
determined from pure noise is almost as large as the values 
retrieved from the dynamo snapshots, with the exception 
of the dynamo run C with a unipolar background field. 
This means that hardly any signal of B^pp remained above 
the noise after spatial smearing and introduction of noise. 
Dynamo C run with a unipolar field is a special case because 
it has a much higher mean field (cf. Figure [2j). The more 
magnetic flux is introduced, the more field can be tangled 
by the turbulent flows, hence noticeably more horizontal 
field is generated. 

The last two rows of Table [3] show the results of the at- 
tempt to estimate how much field the dynamo simulations 
would have to contain to reproduce the observed (B^pp) and 
{\B^pp\) values. Multiplying the magnetic field strengths by 
factors of 2 and 3 (everywhere in the simulation box) in the 
case of the dynamo run G and dynamo run C, respectively, 
gives a mean total magnetic field strength of 170 G and 67 G 
at the levels of r = 1 and r = 0.1, respectively. The strength 
of the mean vertical magnetic field at the same levels be- 
comes 84 G and 27 G, respectively. These values are con- 



sistent with the extrapolation by iPietarila Graham et al.l 
(;2QQ9l )on the basis of Hinode data. Also, as we have seen in 
Fig [2l the average mean field strength distribution of run 
G can be reproduced by multiplying the run C result by a 
factor roughly corresponding to the square root of the ra- 
tio of the total magnetic energies. This suggests that such a 
simple scaling might extend somewhat into the regime 
that is not covered by the simulations. The necessary scal- 
ing factor of 2 for run G is not unreasonably large. 

Figure |4] shows one snapshot from the dynamo run C 
with the field multiplied by a factor of 3, before and after 
spatial smearing, together with an observed region of the 
same size, outlined by the white square in Fig [H Maps of 
the longitudinal B^pp and transversal B^pp magnetic flux 
density are shown. The granular pattern is indicated by 
white contours of normalized continuum intensity equal to 
1.05. Vertical fields with mixed polarities on scales of less 
than a third of an arcsec are smeared into unipolar 'tube' 
or 'sheet '-shaped patches (yellow and blue in the central 
frame). The bundles of horizontal field, composed of loops 
of different sizes, are molded into patches with higher B^pp 
values. These can occur on the edges of granules or be- 
tween patches with vertical field of opposite polarity (see, 
for instance, the feature at [4.5'', 2'']). Similar features can 
be seen in the Hinod e maps, e.g. in t he tw o l ower frames 
on the right (see also llshikawa et al.l (|2QQ8f ) or iLites et aD 
(2008*)). At the position [5'', 3.5''] in the central frame, a ver- 
tical field of B^pp ^ 20 Mx/cm^ is loc ated inside a granule. 
Asimilar case has been observed by lOrozco Suarez etldl 

(goosl). 

Figure [5] shows a different region on the Sun, with more 
vertical flux (outlined by the yellow square in Fig [T]) com- 
pared with a snapshot from the dynamo run C with a 
unipolar background field. No scaling of the dynamo field 
has been carried out here. The maps based on the simula- 
tion show features which are very similar to the structures 
present in the observed weak network region, although the 
B^pp signal is somewhat weaker than in the observations. 
Small patches of horizontal field correspond to small loops 
that are visible between concentrations of the vertical flux 
at original resolution of the simulations. 

5.2. Comparison with the data set II 

IPietarila Graham et al.l (|2QQ9t ) demonstrated how different 
effects can influence the probability density function (PDF) 
of magnetic field derived from Stokes V spectra, such that 
it differs from the PDF of the underlying magnetic field. 
In particular, the effect of noise leads to a PDF which 
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Fig. 4. Comparison of a snapshot from dynamo run C with 
magnetic field scaled by a factor 3 at original (left column) 
and Hinode (middle column) resolution, with Hinode ob- 
servations (right column). From top to bottom: normalized 
intensity, longitudinal and transversal apparent magnetic 
flux density. The observed region is outlined by a white 
square in Fig[Tl White lines outline the contours of normal- 
ized continuum intensity equal to 1.05. 



has a peak at the position that corresponds to the noise 
level. Here we compare the PDFs of B^pp (left-hand panels) 
and B^pp (right-hand panels) determined from the synthe- 
sized Stokes profiles with those from observations. Figure [6] 
shows, from top to bottom, PDFs computed for the mixed 
polarity (20 G) snapshot, a snapshot from the dynamo run 
C without scaling and from the dynamo run C after scaling 
by a factor of 3. Overplotted are the results calculated from 
pure noise with a standard deviation corresponding to the 
noise level of the observations (dotted lines) and the PDFs 
obtained from data set II (red lines). 

The figure shows that all the PDFs from the simulation 
snapshots at Hinode resolution are strongly influenced by 
the noise at the smallest field strengths, which is also the 
case for the PDF based on the Hinode data. The noise- 
induced maxima lie at approximately B^pp ^ 1 G and 
B^pp ^ 20 G. Due to the mixed polarity field on the small 
scales, dynamo run C shows significant loss of stronger 
signals after spatial smearing. The mixed-polarity simula- 
tion snapshot, on the other hand, contains larger unipo- 
lar patches so that it retains a considerable amount of the 
stronger vertical field. Its PDF for B^pp has an extended 
tail of stronger field, which corresponds to the observed 
distribution. In the observations this extended tail is a 
result of the contribution of the network flux concentra- 
tions. The B^pp distributions of both, mixed-polarity and 
dynamo C run snapshots, follow closely the distribution 
generated from the pure noise, which means that the sig- 
nature of the horizontal field is mostly lost in the noise. 
Only a small percentage of the pixels show a B^ signal 



1 2 3 4 5 6 1 2 3 4 5 6 20 21 22 23 24 25 26 
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Fig. 5. Same as Fig. |4]but for the dynamo run C with a 
background unipolar field and a different region from the 
observed map (outhned by a yellow square in Fig [T]) . 



above the noise level. However, after scaling the original 
field values by a factor 3, the PDF for B^pp (bottom row of 
the Fig. [6]) agrees well with the observed one. Small discrep- 
ancies at the smallest signals result from poor sampling (a 
consequence of the small region covered by the simulations). 
The corresponding synthesized distribution of B^pp follows 
the observed distribution up to approximately 10 Mx/cm^. 
Pixels with B^pp signals higher then 40 Mx/cm^ are very 
rare for the maps resulting from dynamo run C. 

5.3. Changes with the heliocentric angle 

The variation of the apparent magnetic flux density at dif- 
ferent heliocentric angles has been determined using a snap- 
shot from the dynamo run C. The magnetic field strength is 
scaled by a factor of 3 in keeping with the findings from the 
previous sections. The upper panels of Fig. [7] show height 
profiles of the mean transversal (right) and mean absolute 
longitudinal (left) field in the simulation at different he- 
liocentric angles. Since the transversal component of the 
magnetic field becomes increasingly dominant in the layers 
above optical depth unity, the mean absolute longitudinal 
(line-of-sight) component of the field increases as the line of 
sight becomes inclined with respect to the surface normal. 
However, the polarization signals in spectral lines reflect the 
component of magnetic field averaged over the line forma- 
tion heights. The lower left panel of Fig. [71 shows the PDF 
of Bavei the vertical component of magnetic field averaged 
over the height range that corresponds to logr = [—3.5, 0.1] 
at different heliocentric angles. The PDFs at fi < 0.8 follow 
closely the PDF at /i = 1, with a discrepancy at the strong 
field end. This discrepancy increases with the heliocentric 
angle. 

As illustrated bv iPietarila Graham et all (|2009f ). Bave 
and \B^ \ are well correlated, which explains the similar 
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Fig. 6. Probability density function (PDFs) for the longitu- 
dinal (left column) and transversal (right column) apparent 
magnetic flux density. PDFs from synthetic Bapp at origi- 
nal (solid) and Hinode (plus signs) resolution are compared 
with the observed PDFs from data set II (red solid line). 
From top to bottom: mixed polarity simulation (20 G), dy- 
namo run C and the same snapshot scaled by a factor 3. 
Also shown are the PDFs derived from Stokes parameters 
resulting from pure white noise with a standard deviation 
of 3 X 10-^ (dotted lines). 



trend of ^^^^ heliocentric angle, as shown in the 

lower right panel of Fig. [71 Plotted are the mean values of 
\B^pp\ and B^^^ as functions of the heliocentric angle. The 
mean value of \B^pp\ stays roughly constant until /i = 0.8 
and then gradually decreases when moving further away 
from disk center. The mean value of decreases mono- 
tonically. 

5.4. Comparison with observations at /j. = OA (6 = 66"^ ) 

In this section, we compare the synthesized signals at 
Hinode resolution with the observations. The upper panel 
of Fig. [8] shows PDFs of the field proxies, B^pp and 5^^, 
retrieved from the observations at disc center and 6> = 66° 
(data sets I and III), respectively. The maxima of the PDFs 
are displaced with respect to each other, which is due to 
the difference in the noise level (Table [T]) as indicated by 
the overplotted PDFs resulting from the pure noise. The 
probability of occurrence of the B^pp signals larger than 
30 Mx/cm^ decreases at 6> = 66°, while the probability 
of occurrence of the B^pp signals larger than 200 Mx/cm^ 
increases. This effect results from the strong magnetic con- 
centrations present in the field of view. When a strong mag- 
netic tube-like structure with a vertical field of kG strength 
is observed under an angle then B^pp signal decreases, while 
B^pp signal displays a corresponding increase relative to the 
value at the disc center. This is because, projected on to 
the line of sight, the strong vertical component of magnetic 
field gives rise to a significant transversal component, lead- 
ing to a strong B^ signal. Note, however that more pixels 




B^pp [Mx/cm^] B^„,, [Mx/cm^] 



Fig. 8. PDFs for the longitudinal (left column) and 
transversal (right column) apparent magnetic fiux density. 
Upper row: Bapp distributions obtained from observations 
at heliocentric angles ^ = 0° (red solid line) and = 66° 
(green dashed line). The dotted lines in the same colors give 
the PDFs resulting from the corresponding noise. Bottom 
row: Comparison of observations with the synthetic Bapp 
from dynamo run C (scaled with a factor 3) at heliocentric 
angles 6 = 0° (black plus signs) and 6 = 66° (blue plus 
signs), at Hinode resolution. 

has values in the range 70 < B^pp < 200 Mx/cm^ where 
FDF{B^pp) decreases towards the limb. Hence the mean 
value of B^pp decreases with heliocentric angle, in agree- 
ment with the result obtained from local dynamo simula- 
tions at original resolution (shown in lower right frame of 
Fig.©. 

The lower panels of Fig.[8]show the PDFs retrieved from 
the dynamo run C snapshot at Hinode resolution. Their 
trends correspond to the ones at original resolution, show- 
ing a small difference only for the strong signals. The PDF 
of the synthesized B^pp signal follows the observations at 
both heliocentric angles up to B^pp ^ 200 Mx/cm^. The 
synthesized B^pp PDF matches the observations at 6> = 66° 
up to approximately 20 Mx/cm^. 

6. Discussion 

We have based our comparison between MHD simula- 
tions and observational results upon the magnetic prox- 
ies, the longitudinal and transversal appar ent magnetic 
fiux density introduced bv iLites et al.l (|2008f ). It has been 
shown that these proxies are prone to the infiuence of the 
non-magnetic properties of the atmosphere (|Beck fc Rezaeil 
|2009[ ) and do no t refiect correctly the proper ties of the 
underlying field (|Pietarila Graham et al.l [2009f ). However, 
if we take that realistic MHD simulations represent the 
solar atmosphere properly (supporte d by nu merous stud- 
ies, e.g.iS chiissler et al. (2003); ShelyaJetaD (2004, 2003); 
iDanilovic et al.l (^2009)) and that our forward modelling of 
the instrumental effects and the noise is correct, we intro- 
duce the same assumption as in the case of the observations, 
allowing us to make a rela t ively unbiased comparison with 
the findings of iLites et all (|2008D . 

The magnetic proxies obtained from the synthesized 
data are compared with observations of the internetwork 
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Fig. 7. Simulated change with heliocentric angle retrieved from a dynamo run C snapshot. Upper panels: The mean 
absolute longitudinal (left) and transversal (right) components of the magnetic field as a function of height in the 
simulation box at different heliocentric angles. Vertical lines indicate the level of optical depth unity. Lower left: PDFs 
for the signed longitudinal field (averaged over the line formation height) at different heliocentric angles. Lower right: 
Change of the mean absolute longitudinal (l^^ppl) mean transversal (5^^) apparent flux density as a function of 
the heliocentric angle. These values are calculated from synthesized Stokes profiles. 



region. The mixed polarity simulations (which have a mag- 
netic Reynolds number below the threshold for dynamo ac- 
tion) reproduce the observed v ertical flux den sity, as was 
previously shown by iKhomenk o et al.l (|2QQ5a[ ). However, 
they do not contain enough horizontal fleld, to be consistent 
with the Hinode data. Simulations of surface dynamo give 
a ratio of the horizontal to the vertic al flux density con- 
sisten t with the observational results (|Schussler fc VogTerl 
l2QQ8f ), but the overall amplitudes are too low, at least for 
the simulations with the magnetic Reynolds number feasi- 
ble at the moment. A simple scaling of the fleld in the sim- 
ulation domain by a flxed factor brings the magnetic flux 
density from the dynamo snapshots at Hinode resolution 
close to the observed ones. Some justiflcation for this ten- 
tative procedure comes from comparing simulations with 
different magnetic Reynolds numbers. The resulting mean 
magnetic fleld strength is of the order of 70 G at r = 0.1 
and around 170 G at the solar surf ace. This is roughly 
consis tent with the estimates given bv lTruiillo Bueno et al.l 
(120041) . since the line they consider samples the middle pho- 
tosphere, a nd with extrapolations based on the cancellation 
coefflcient ([Pietarila Graham et al.n2009i) . 



The contribution of the surface dynamo dominates 

< 30 Mx/cm^ and < 



the 



regions 



with \BtJ 



200 Mx/cm at Hinode spatial resolution. Everything 
stronger than th at could imply: (1) fl ux emergence with kG 
horizontal flelds (|Cheung et al.|[2QQ8l ) or (2) strong network 
flelds, which have a different source from those in the in- 
ternetwork, a possibility supported by studies of ep hemeral 
active regions (|Harvevlll993l : iHagenaar et al.lll999[ ). 

The synthesized magnetic proxies exhibit a decrease 
from the disk center to the limb. This is in quali tative 
agreement with the observations bv lLites et al.l (|2008f ). The 
PDFs show little /i dependence for the small values of 
the magnetic proxies, in accordance with t he results of 
iMa rtinez Gonzalez et al.l (|2008t ). In contrast, ISteiner et al.l 
(20^ could not reproduce the observed simultaneous de- 
crease towards the limb of the longitudinal and transversal 
fleld proxies on the basis of MHD simulations assuming a 
prescribed vertical magnetic flux or a transport of horizon- 
tal flux into the computational box from below. 

Our results support the presence of local dynamo action 
in the quiet Sun. They also suggest that the magnetic fleld 
is dominantly horizontal. The ratio of the mean horizontal 
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and vertical component of the dynamo-generated magnetic 
field reaches values between 2 and 4 in the optical depth 
interval — 2 < logr < —1. This gives the observed ratio 
of the transversal and longitudinal apparent magnetic fiux 
densities at the Hinode resolution when an adequate noise 
level is considered. 

We note, however, that the conclusions drawn here are 
based on a simple assumption drawn from the properties 
of magnetic field in the simulations with different magnetic 
Reynolds numbers. The increase in the field strength by a 
factor of 2 or 3 would result in a change of physical pa- 
rameters as well as the dynamics of the plasma since the 
strong field is much less susceptible to the infiuence of the 
turbulent motions. These effects are not taken into account 
in our scaling procedure. Thus, we consider our work to be 
a step towards revealing the role of the surface dynamo in 
the quiet Sun magnetism until more realistic simulations, 
as well as higher resolution observation of the solar photo- 
sphere become available. 

Acknowledgements. We thank Jonathan Pietarila Graham for the 
valuable help with the PDFs and Anna Pietarila for a critical reading 
of the paper. Hinode is a Japanese mission developed and launched by 
ISAS/JAXA, with NAOJ as domestic partner and NASA and STFC 
(UK) as international partners. It is operated by these agencies in 
co-operation with ESA and NSC (Norway). This work was partly 
supported by WCU grant No. R31-100016 from Korean Ministry of 
Education, Science and Technology. 



Orozco Suarez, D., Bellot Rubio, L. R., del Toro Iniesta, J. C, & 

Tsuneta, S. 2008, A&A, 481, L33 
Pietarila Graham, J., Danilovic, S., & Schiissler, M. 2009, ApJ, 693, 

1728 

Petrovay, K., & Szakaly, G. 1993, A&A, 274, 543 

Sanchez Almeida, J., Emonet, T., & Cattaneo, F. 2003, ApJ, 585, 536 
Schiissler, M., Shelyag, S., Berdyugina, S., et al. 2003, ApJ, 597, L173 
Shelyag, S., Schiissler, M., Solanki, S. K., Vogler, A. 2007, A&A, 469, 
731 

Shelyag, S., Schiissler, M., Solanki, S. K. et al., 2004, A&A, 427, 335 
Solanki, S. K. 2009, Astronomical Society of the Pacific Conference 
Series, 405, 135 

Steiner, O., Rezaei, R., Schaffenberger, W., &; Wedemeyer-Bohm, S. 

2008, ApJ, 680, L85 
Steiner, O., Rezaei, R., Schlichenmaier, R., et al. 2009, 

arXiv:0904.2030 

Piskunov, N. E., Kupka, F., Ryabchikova, T. A., et al. 1995, A&AS, 
112, 525 

Thevenin, F. 1989, A&AS, 77, 137 

Trujillo Bueno, J., Shchukina, N., &: Asensio Ramos, A. 2004, Nature, 
430, 326 

Tsuneta, S., Ichimoto, K., Katsukawa, Y. et al. 2008, Sol. Phys., 249, 
167 

Vogler, A., Shelyag, S., Schiissler, M., et al. 2005, A&A, 429, 335 
Vogler, A. 2003, PhD Thesis, University of Gottingen, Germany, 

http: / / webdoc.sub.gwdg.de / diss / 2004 / voegler 
Vogler, A., & Schiissler, M. 2007, A&A, 465, L43 
Schiissler, M., & Vogler, A. 2008, A&A, 481, L5 

de Wijn, A. G., Stenfio, J. O. Solanki, S. K. and Tsuneta, S. 2008, 
Space Sci. Rev. , Volume 144, Issue 1-4, pp. 275-315 



References 

Asensio Ramos, A. 2009, ApJ, 701, 1032 

Asensio Ramos, A., Martinez Gonzalez, M. J., Lopez Ariste, A. et al. 

ApJ, 659, 829 
Beck, C, & Rezaei, R. 2009, A&A, 502, 969 

Bello Gonzalez, N., Yelles Chaouche, L., Okunev, O., &; Kneer, F. 

2009, A&A, 494, 1091 
Berdyugina, S. V., & Fluri, D. M. 2004, A&A, 417, 775 
Cattaneo, F. 1999, ApJ, 515, L39 
Centeno, R., et al. 2007, ApJ, 666, L137 

Cheung, M. C. M., Schiissler, M., Tarbell, T. D., &; Title, A. M. 2008, 
ApJ, 687, 1373 

Danilovic, S., Gandorfer, A., Lagg, A., et al. 2008, AfcA, 484, L17 
Danilovic, S., Schiissler, M., &; Solanki, S. K. 2009, arXiv:0910.1211J 
Dommguez Cerdefia, I., Sanchez Almeida, J., Kneer, F. 2003, A&A, 

407, 741 

Frutiger, C, Solanki, S. K., Fligge, M., & Bruls, J. H. M. J. 2000, 
A&A, 358, 1109 

Hagenaar, H. J., Schrijver, C. J., Title, A. M., & Shine, R. A. 1999, 

ApJ, 511, 932 
Harvey, K. L. 1993, Ph.D. Thesis, Utrecht University 
Harvey, J. W., Branston, D., Henney, C. J., & Keller, C. U. 2007, 

ApJ, 659, L177 
Ishikawa, R., et al. 2008, A&A, 481, L25 
Kosugi, T., et al. 2007, Sol. Phys., 243, 3 

Khomenko, E. V., Collados, M., Solanki, S. K., et al. 2003, A&A, 

408, 1115 

Khomenko, E. V., Martinez Gonzalez, M. J., Collados, M., et al. 2005, 
A&A, 436, L27 

Khomenko, E. V., Shelyag, S., Solanki, S. K., Vogler, A. 2005, A&A, 
442, 1059 

Lites, B. W., Elmore, D. F., & Streander, K. V. 2001, in: ASP conf. 

ser.. Vol 236, Advanced Solar Polarimetry - Theory, Observation, 

and Instrumentation, ed. M. Sigwarth, 33 
Lites, B. W., et al. 2008, ApJ, 672, 1237 
Martin, S. F. 1988, Sol. Phys., 117, 243 

Martinez Gonzalez, M. J., Collados, M., Ruiz Cobo, B., & Solanki, 

S. K. 2007, A&A, 469, L39 
Martmez Gonzalez, M. J., Asensio Ramos, A., Lopez Ariste, A., &; 

Manso Sainz, R. 2008, A&A, 479, 229 
Orozco Suarez, D., et al. 2007, ApJ, 670, L61 
Orozco Suarez, D., et al. 2007, PAS J, 59, 837 



